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Abstract 

Using the self-consistent Ornstein-Zernike approximation (SCOZA) results 
for the 3D Ising model, we obtain phase diagrams for binary mixtures described 
by decorated models. We obtain the plait point, binodals, and closed-loop coex- 
istence curves for the models proposed by Widom, Clark, Neece, and Wheeler. 
The results are in good agreement with series expansions and experiments. 
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1. Introduction 



Critical phenomena of mixtures are considerably more complicated to inves- 
tigate than such phenomena in pure fluids. Puriflcation of components, precise de- 
termination of composition, gravitational effects and long equilibration times make 
experiments on mixtures much more difficult to realize. Similar difficulties are found 
in theoretical analyses, where the addition of one component greatly complicates the 
generalization of the lattice-gas model to mixtures, for example the derivation and 
analysis of series expansions 

The use of decorated lattice-gas models to describe fluid mixtures represents 
an important gain in precision and information about these models. This is a result 
of the exact mapping between these models and the single-species lattice-gas or Ising 
model. The exact results in two dimensions, and exhaustive series analysis in three, 
allow a detailed description of the phase diagrams, coexistence surfaces and critical 
exponents of the much more complicated models for mixtures. 

In this paper we obtain information related to critical phenomena of binary 
mixtures for various decorated lattice models using Ising-model results obtained via 
the SCOZA [0]. This approximation provides very precise results for the lattice-gas 
Ising model in three dimensions as shown elsewhere 0. 

The next sections deflne the models we work with and show how to derive the 
grand partition function for these models. Then we present a brief discussion about 
SCOZA. Finally, we present the phase diagrams and coexistence curves we obtain for 
the binary mixture models. 

2. Decorated model 

Consider a binary mixture of A^i molecules of type 1 and N2 molecules of type 
2 contained in a volume V . A decorated lattice model for a mixture is constructed as 
follows We divide the volume into C primary cells, which are centered on sites 

of a regular lattice, and secondary cells centered on the bond midpoints of this lattice. 
If the primary cells have volume fo, then the secondary cells occupy a volume 2fo/g, 
where q is the lattice coordination number. Each primary cell overlaps q secondary 
cells and each secondary cell overlaps two primary cells, as shown in flgure 1. A 
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primary cell is occupied if a type-1 molecule lies within its boundaries; secondary cells 
are occupied by type-2 molecules. The cells are not meant to restrict the molecules 
to lattice sites. They merely provide a discrete coordinate system for measuring the 
positions of particles of type 1 and 2 and thus defining intermolecular separations. 

[Insert figure 1 about here] 

All interactions are pairwise. If two molecules of type 1 are in the same primary 
cell, the interaction energy is +oo; this energy is — e if they occupy adjacent primary 
cells, and zero otherwise. The interaction energy for a pair of type-2 molecules is 
+00 if they occupy the same secondary cell and zero otherwise. If a primary cell 
containing a type-1 molecule overlaps a secondary cell occupied by a type-2 molecule, 
the interaction energy is 0. Particles of different species have no interaction if their 
cells do not overlap. 

The condition that there is no interaction between type-2 molecules in different 
secondary cells is essential to the existence of an exact mapping between the decorated 
model and the spin-| Ising model. As a consequence, the pure component 2 cannot 
undergo a phase transition 

So far we have described the decorated model in general terms. The specific 
models for binary mixtures by Widom, Clark, Neece, and Wheeler are obtained by 
adjusting the parameters for the interaction energy. 

Widom [^] defined the simplest model by setting a hard-core repulsion be- 
tween molecules of types 1 and 2, the only interaction beyond exclusion of multiple 
occupancy. The model undergoes a first-order phase transition even though the in- 
termolecular forces are all assumed to be infinitely strong repulsions. 

Clark generalized Widom's model by allowing a finite repulsion between 
unlike nearest-neighbour molecules. All other interactions between like molecules are 
zero and multiple occupancy of cells is not allowed. This model has a coexistence 
dome with a maximum critical solution temperature. 

Neece ^ extended Widom's model including the effects of interactions between 
molecules on the primary lattice. This model exhibits a variety of behaviours with 
changing temperature, especially if the molecules on the primary lattice interact with 
an attractive potential. It provides a good qualitative description of certain types of 
'gas-gas equilibrium' such as that found in the He-Xe system. 
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Wheeler ||^ studied the binary case with an attractive interaction between 
molecules on the primary lattice, and a finite interaction (attractive or repulsive) 
between type-1 and -2 molecules. This version proved useful in understanding critical 
phenomena in very dilute solutions. 

Another model studied by Wheeler ||^ considers a highly directional, short- 
ranged interaction favoring special orientations of the molecules at low temperatures. 
Unlike the other models, a primary or secondary lattice site can be occupied either by 
a type-1 or a type-2 molecule. Interactions between particles of the same species (1-1 
or 2-2) are zero regardless of their orientations. The 1-2 interaction depends upon the 
orientation of the particle in the secondary cell. If the 1-2 bond involves the special 
contact point of this particle then the energy is negative (attractive) otherwise it is 
positive (repulsive). Each molecule (1 or 2) has one special contact point and {u — 1) 
other identical contact points regardless of whether it is on a primary or secondary 
cell. We consider the case with uj = Q. Thus, if a molecule of either type is in a 
secondary cell which has one l-fiUed and one 2-filled primary neighbour, then it has 
one orientation with energy U2 < and five orientations with energy Ui > 0. If both 
their primary neighbours are of the opposite type, then there are two orientations in 
which it makes one attractive and one repulsive contact with resulting energy U2 + U1, 
and four orientations in which it makes two repulsive contacts with resulting energy 
2Ui. The 1-2 interaction energy is independent of the orientation of the particle in the 
primary cell. This is required to map the model exactly onto the spin-1/2 Ising model. 
This model presents closed-loop coexistence curves with both upper and lower critical 
solution temperatures, analogous to the behaviour found in the nicotine + water and 
m-toluidine + glycerol systems. 

Table 1 shows a summary of the models presented above. 

[Insert table 1 about here] 

3. Derivation of the Grand Partition Function 

In this section we derive the grand partition function for the decorated models 
and show how to map these models into the lattice gas Ising model. We follow the 
derivation given by Wheeler (for details see |]I|). The grand partition function can be 
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written as 

Z = J2I1 z^'x^'^'z^'e-^ (1) 

Ci Ca 

where the sums are over all assignments of filled or empty primary and secondary 
cells. is the sum of all primary-secondary interactions, A^n is the number of 

fiUed-fiUed nearest neighbour pairs of primary cells and x — exp(— e/ZcT). zi and Z2 
are dimensionless activities defined as, 

Zl = voi^ J e-, (2) 

fqvo\ / 27rm2kT\ ^ tz 

where d is the dimensionality. 

Let rjij — e"'^*^/*:^^ i^j — 0, 1, where ■000, V'oi, V'li ^^e the energies for a molecule 
of type 2 that lies between two empty, one empty and one filled, or two filled primary 
cells, respectively. The idea is to sum over filled and empty states for each secondary 
cell, once a specific configuration of the primary cells is chosen. Thus if a secondary 
cell lies between empty-empty primary cell pairs, it contributes to the grand parti- 
tion function a factor 1 + -22^001 if it li^s between empty-filled primary cell pairs, it 
contributes a factor 1 + Z21]qi, and if it lies between filled- filled primary cell pairs, it 
contributes a factor 1 -|- Z2r)u- Then, 

Z = + ^277oo)^°°(l + Z2V0lf°'{l + Z2Vuf'' (4) 

Ci 

where A^oo and A^oi are the number of empty-empty and empty-filled nearest neighbour 
primary cells. Using the lattice identities, qNi — 2Nii + Nqi, qNq — 2A^oo + -^01, -^1 = 
N and Nq ^ C - N we get 

Z = (1 + ^2?7oo)^ l^^^exp (^^) = (1 + ^2^00) ^^i?- (5) 
Comparing this with the partition function for the lattice gas, 

2r = J2 {zufexpiPneRNu), 
c 
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we obtain the relations that map the decorated model exactly to the lattice gas Ising 
model: 



qC 

Z = {1 + Z2r]m) 2 Zr 



Zr = Zi{l + Z2ri00) (l + 22^?0l)' 

exp(-^) = x(l + Z2'7oo)(l + 2;2?7ii)(l + 22'7oi)~^- (6) 



^kTRy 

The subscript R denotes that the quantity refers to the reference system. This is 
an exact correspondence between the decorated model and the spin- 1/2 Ising model. 
All thermodynamic properties can be obtained from the grand canonical potential 
'n'{zi,Z2,x) = limc_-»oo 1/Cln^ by differentiation. 

The thermodynamic pressure p{fi, T) is given by 

1 q 
pp = lim — In Z{T, ij,,V) = - In (1 + Z2'noo) + Pr Pr (7) 

for the decorated model. The number densities are given by: 

P2 = (j^) =U + qiC-A)pR + iA + B-2C)^ (8) 

where $h is the energy density in the reference system and, 

A= ^ , B= ^ , C= , (9) 

(1 + Z2r]oo) (1 + Z2'nu) (1 + Z2r]oi) 

Other quantities of interest are the internal energy, given hj U = {dv: / d[3)zi,z2^ and 
the specific heat at constant volume, given by Cv,x = {dU/dT)pj^^p^. 

4. Self-Consistent Ornstein-Zernike Approximation(SCOZA) 

The SCOZA closure of the Ornstein-Zernike equation yields critical parame- 
ters of remarkable accuracy for three-dimensional lattice gases with nearest-neighbour 
interactions p|, ||. 

In the SCOZA, the direct correlation function c(r) is determined in a way that 
results in a unique free energy function. The starting point is the Ornstein-Zernike 
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relation 

hi = Ci + p ^ Ci_jhj (10) 

Lattice 

where h is the total correlation function. Thermodynamic consistency is embodied in 
the exact relation that must be satisfied by the derivation of the free energy 

d\pu) _ d\(3p) 
^ 9p2 d(3dp ^ ' 

where u is the internal energy per particle, p is the pressure and p the density. Cq = c(0) 
is fixed by the core condition on the total correlation function, ho = —1. The single 
approximation of the method is the truncation of the direct correlation function c(r) 
at the nearest-neighbour separation: c(r) = for r > 1. From the relations between 
Co,Ci and hi we find a partial differential equation for ci(/9, p) that can be solved 
numerically. For details see references 0, ^. 

5. Results 

In this section we specify the mapping equations for each model of interest 
and we also show the corresponding phase diagrams. We study binary mixtures in 
three dimensions. Thus the reference one-component lattice gas has a range of ther- 
modynamic states corresponding to liquid- vapor coexistence terminating in a critical 
point. At the boundary of the two-phase region in the N/Q) plane of the refer- 

ence system the activity zr, the temperature Tr, and the pressure-temperature ratio 
Pr/Tr are all singular, having discontinuous derivatives there. Therefore, we see from 
equations ( ^ and (|^), that the functions zi,Z2,p/T of the two-component models are 
singular at the corresponding point in the {pi, P2) plane. 

As a reminder, in our notation e is the interaction energy between two type 1 
molecules when sitting at neighbouring primary cells, and (p denotes the interaction 
energy between a type-1 and a type-2 molecule when they occupy overlapping primary 
and secondary cells. 

5.1. Reference system 

Our reference system is the lattice gas transcription of the Ising model, which 
has two-phase behaviour below a critical temperature Be- The coexistence curve is 
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shown in figure 2. In the two-phase region there are two values of the density which 
are in equihbrium for each temperature 9. The same is true for the binary system 
after carrying out the transformations mapping it to the reference system. For each 
temperature below the critical temperature in the binary system, the binodal curve 
corresponds to the coexistence curve in the reference system. By varying the temper- 
ature the binodals will then form the coexistence surface. 

[Insert figure 2 about here] 

Before investigating the nature of the coexistence surface it is important to re- 
view some of the properties of the one-component lattice gas in the two-phase region. 
Nearest neighbour lattice gases of dimensionality greater than one have a critical den- 
sity pc = 1/2) with the critical isochore corresponding to this density. The properties 
along the critical isochore are known exactly only in the two-dimensional case. How- 
ever, the behaviour in three dimensional systems is known with sufficient accuracy for 
our purposes. In this paper we use the results obtained via SCOZA. 

The equation of the coexistence curve can be written as ^ = t(p), for which 
there are two values of p for each 9 below the critical temperature 9^- The activity 
along the critical isochore is 

A = exp . (12) 

As 9 ^ 9c in the two-phase region, the usual power law expressions are assumed 
to hold. The coexistence curve in the neighbourhood of the critical density can be 
written as 

9c-rip)r^\p-pc\^, (13) 

where the exponent 1//? is very close to 3 in three dimensions, so the coexistence curve 
is nearly cubic. In SCOZA, l/f] = 20/7 [Q. The heat capacity (constant volume) 
diverges as the critical point is approached asymptotically along the critical isochore, 
following 

r~(^^,-e)-"', (14) 

where a' > and nearly equal to zero in three dimensions. In SCOZA, a' = 1/10 
|]. Current best estimates for (3 and a' are 0.3257 ± 0.0025 and 0.1130 ± 0.0075 ig, 
respectively. 
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5.2. Behaviour of a general binary model 



The complete three-dimensional phase diagram for a binary system at constant 
pressure is shown in figure 3. 

[Insert figure 3 about here] 
The coexistence surface separates the one- and two-phase regions of the system. 
In a constant temperature plane the variables are the number densities of the two 
components. The two-phase region is bounded by a binodal curve APB, containing 
tie lines between points representing the compositions of the two coexisting phases. 
The tie lines vanish at point P, called the plait point (isothermal critical mixing 
point), where the coexistent phases lose their separate identities to merge into a 
single homogeneous phase. In the full three-dimensional diagram the plait points 
form a curve on the coexistence surface. This curve ends up reaching a maximum 
at the top C, above which only one phase is found. C is then a critical point of the 
system. We will see that in some cases the binodals are closed loops involving two 
plait points. 

5.3. Widom's model 

The number densities for Widom's model, which is defined by e = and 
= -l-cxD, are 

Pi = PR 

P2 = 77 ( ~ + ~ ) [exp(-/3ijeif) - 1]. (15) 



3 V2 tR 

In order to obtain the binodal curve, shown in figure 4, we use the reference sys- 
tem coexistence curve values for the energy density, temperature and number density 



in equation 15 



[Insert figure 4 about here] 
The binodal curve APB and the line AB for T = separates the two-phase 
region from the remainder of the composition triangle corresponding to the homoge- 
neous states of the system. In order to determine the tie line in the binary system we 
identify the double-valued density in the reference system for a given temperature. 
The corresponding densities in the binodal are coexistent. The tie line CD connects 
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pairs of points on the binodal curve which represent the compositions of coexistent 
phases. The longest tie hne connects the coexistent densities (pi = 0,p2 = 1) and 
(pi = l,p2 = 0), corresponding to 6' = in the reference system. The tie hnes vanish 
at the point P, where the conjugate phases merge into a single homogeneous phase. 
Thus P is a critical point achieved by changing composition at fixed temperature. P 
is frequently called a plait point in order to distinguish it from the usual critical point 
achieved by changing temperature. The coordinates Pi,P2 of the plait point are 

pt = i p; = ^ (^) [exp{-(3,en) - 1] = 0.1968242 (16) 



The binodal curve is an isothermal section of a more general three-dimensional 
diagram. For this model the number density and all the other quantities derived are 
independent of the temperature T. This happens because the molecular interaction 
potential is either or +oo, so the total potential energy of this system never deviates 
from 0, i.e. there is no energy scale. Nevertheless, the model exhibits phase coexistence 
and an associated plait point. 

5.4. Clark's model 
For this model, e = and (j) is finite. The number densities pi and p2 are 



Pi = PR 
1 

P2 = - 



|a + iv|(C -A) + Nn{A + B- 2C) 



(17) 



where A, B and C were defined by equation ^ when deriving the grand partition 
function. 

Equations ^ and |l^ are the mapping between Clark's model and the lattice gas 



reference system. In order to obtain the binodal curve, we determine the activities 
and the number densities for the binary model by using the energy density, number 
density and temperature values of the reference system coexistence curve. 

[Insert figure 5 about here] 
Figure 5 shows the isothermal binodal curves in the (pi, P2) plane, correspond- 
ing to (j) = 0.5 and temperatures T = 0.1; 0.2; 0.247162. The coexistence surface, 
figure 3, is the set of all such binodal curves, which close to form a dome whose apex 
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marks the critical point (Tc = 0.247173). Each isotherm in the {pi, P2) plane, has an 
upper and a lower plait point. These plait points form a plait point curve that ends 
at the critical point Tc. 

5.5. Neece's model 
For this model e is nonzero and = +00. The number densities are 
Pi = PR 



P2 = ^(l + iVg + iVn) 



1 , e/J e 



If e < there is an attraction between nearest neighbours on the primary lattice. In 
our model \tB\ = \e\ = 1. The second term in brackets in the expression for p2 is never 
zero for T > Be- If T = 6c, then this factor is zero and p2 = at pi = 1/2. If T < 6c, 
the second factor is zero for some value of 6 which in turn corresponds to two values 
of pi. The binodal curve for e < is shown in figure 6 at different temperatures. 
The critical temperature is Tc = 1.1299, represented by curve C. From the inset of 
figure 6, which is a zoom of curve D, we can see that there are two separate regions 
corresponding to the homogeneous states of the system for T < 6c. Also, in this case 
there is no plait point since there is no tie line of vanishing length. 

[Insert figure 6 about here] 
For e > 0, the binodals as a function of temperature are shown in figure 7. The 
plait points lie along pi = 1/2. 

[Insert figure 7 about here] 

5.6. Wheeler's model 

For this model, e < and (p is finite (attractive or repulsive). The number 
densities pi and p2 are given by 



Pi = PR 

P2 = {hA + q{C-A)pn + {A + B-2C)^ (19) 



where A, B and C are defined by equation |^. 
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[Insert figure 8 about here] 
Figure 8 shows the binodal curves for $ < and $ > at different tempera- 
tures. Curve A is for T = 1.11, curve B is for T = 1.25 and curve C for T = 1.336 for 
attractive and repulsive interactions. The slope of the tie lines is positive for $ < 
and negative for $ > 0. At T = Tc = 1.336498 the binodal shrinks to a point. In a 
3-D diagram, including the temperature, it corresponds to a dome that separates the 
two-phase region from the homogeneous one. 

5.7. Wheeler's oriented model 

In this model particles of the same species (1-1 or 2-2) do not interact. The 1-2 
interaction depends upon the orientation of the particle in the secondary cell. If the 
1-2 bond involves the special contact point of this particle then the energy is [/2 < 
(attractive), otherwise it is t/i > (repulsive). 

[Insert figure 9 about here] 

Figure 9 shows the coexistence curve for r = \U2\/Ui = 0.189 (curve A), 0.18 
(curve B), and 0.17 (curve C). Each curve forms a closed loop with a lower and upper 
critical temperature. The ratios Tu/Ti of upper to lower critical solution temperature 
are respectively 1.161, 1.467 and 1.721 for these curves. Our results are in close 
agreement with those of Wheeler [|rT|, using low temperature series expansions and 



with the experimental data. 

In figure 10 we show the experimental coexistence curves for solutions of glyc- 
erol and guaiacol (circles) [0 and m-toluidine (squares) [jl3[. These curves are qual- 
itatively similar to the closed loops obtained from Wheeler's oriented model. 

[Insert figure 10 about here] 

6. Concluding remarks 

We study binary systems through various versions of the decorated lattice 
model using the self-consistent Ornstein-Zernike approximation. We obtain the bin- 
odals, plait points, and the critical temperatures for the models introduced by Widom, 
Neece, Clark and Wheeler. SCOZA provides a very useful way of getting accurate 
results for binary systems that can be mapped to the lattice-gas Ising model. 
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We also study Wheeler's oriented model, obtaining its coexistence curve for 
various interaction energies. The results are consistent with experimental data and 
low temperature series expansions. 

In order to analyse the algebraic degree of the binodal at the plait point we can 
use the idea proposed by Widom [Q. It consists in measuring the vertical distance 
Ap2 from the binodal curve to a line drawn tangent to the binodal at the plait point. 
By combining the number densities with the energy density and the coordinates of 
the plait point, we find that 

1 jl-a') 

Ap2 ~ |pi - -1^^. (20) 

The modification of the critical exponent by the factor (1 — a') is an example of 
'renormalization' discussed by Fisher |T^. The SCOZA values for /3 and a' are 7/20 
and 1/10, respectively, providing 18/7 for the algebraic degree of the binodal. 

These results can be extended to study ternary systems as well as pure fluids. 
It is not difficult to derive the energy, heat capacity, pressure for all the models studied 
in this paper. 
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Table 1: Interactions between molecules in adjacent cells for the different decorated 
models. 



Model 


Energy between type 1 molecules 


Energy between unlike molecules 


Widom 


e = 


= +0O 


Clark 


e = 


> or < 


Neece 


e > or e < 


= +00 


Wheeler 


e > 


> or < 


Wheeler oriented 


e = 


C/i > or C/2 < 



Figure Captions 

Figure 1. Example of a decorated lattice in two dimensions. The solid squares rep- 
resent the primary cells while the dashed diamonds represent the secondary cells. 

Figure 2. Coexistence curve for the lattice gas Ising model using SCOZA results. 

Figure 3. Coexistence surface for a binary mixture using the results for Clark's model. 
The top C marks a critical point. 

Figure 4. Binodal curve for Widom's model. CD is a tie line and P is the plait point. 
The composition variables pi and p2 are the number densities of the two components. 

Figure 5. Binodal curve for Clark's model for attractive ($ > 0) and repulsive ($ < 
0) interactions. The outer binodal corresponds to T = 0.1; the intermediate one 
corresponds to T = 0.2 and the inner curve is for T = 0.247, very close to Tc = 
0.247173. 

Figure 6. Binodal curve for Neece's model with e < 0. Curve A corresponds to 
T ^ oo; curve B is for T = 2.0 > T^; the critical temperature Tc = 1.12989 is 
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represented by curve C; and T = 1.0 < Tc at curve D. The inset is a zoom of curve 
D. 



Figure 7. Binodal curve for Neece's model with e > 0. Curve A corresponds to T = 0; 
curve 5 is for T > 0; T — > oo is represented by curve C. 

Figure 8. Binodal curve for Wheeler's model for repulsive < 0) and attractive 
> 0) interactions. Curve A is for T — 1.11, curve B is for T — 1.25 and curve C 
for T — 1.336 for both diagrams. The critical temperature is T = Tc = 1.336498. 

Figure 9. Coexistence curves for Wheeler's oriented model in the temperature-mole 
fraction diagram. Curves A, B and C correspond to r = 0.189, 0.18, 0.17, respectively. 

Figure 10. Temperature-mass fraction coexistence curves for solutions of glycerol with 
guaiacol (circles) and m-toluidine (squares). Temperature is in Celsius. 
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